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An approach for evaluation of observables in analytic versions of QCD 
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We present two variants of an approach for evaluation of observables in analytic QCD models. 
The approach is motivated by the skeleton expansion in a certain class of schemes. We then evaluate 
the Adler function at low energies in one variant of this approach, in various analytic QCD models 
for the coupling parameter, and compare with perturbative QCD predictions and the experimental 
results. We introduce two analytic QCD models for the coupling parameter which reproduce the 
measured value of the semihadronic r decay ratio. Further, we evaluate the Bjorken polarized sum 
rule at low energies in both variants of the evaluation approach, using for the coupling parameter the 
analytic QCD model of Shirkov and Solovtsov, and compare with values obtained by the evaluation 
approach of Milton et al. and Shirkov. 

PACS numbers: 12.38.Cy, 12.38.Aw,12.40.Vv 



Consider an observable 0(Q 2 ) depending on a single space-like scale Q 2 (= 
expansion for this observable exists: 



> and assume that the skeleton 



r°° rit °° " r 00 At 

OsUQ 2 ) = / ^F^(i)a pt (te c Q 2 ) + ^ S ^ 1 n/ T^Ptfe^Q 

J 1 n=2 , = 1 J h 



Fq^i,. . .,t„). 



(1) 



The observable is normalized such that 0{Q 2 ) — a p t at first order in perturbation theory. The characteristic functions 
Fq are symmetric functions and have the following normalization: 



dt 



F$(t) = 1, 



dh dt 2 ^ A 



F£(ti,t2) = 1, 



(2) 



and sY are the skeleton coefficients. The perturbative running coupling a pt (Q ) = a(Q )/n obeys the renormalization 
group (RG) equation: 



da pt (Q 2 ) 
dlogQ 2 



PiaUQ 2 ) 



(3) 



X 



In QCD, the first two coefficients (3q = (1/4)(11 — 2n//3) and fix = (1/16)(102 — 38n//3) are scheme-independent 
in mass-independent schemes; rif is the number of active quarks flavors. The value of C depends on the value of 



the scale A in a 



pt 



A 2 p' 

yi (o) e 



') I n MS scheme C = C = —5/3. The skeleton integrands and integrals are 



independent of C. Expansion Q exists in QED if one excludes light-by-light subdiagrams HQ. In the QCD case, 
the leading skeleton part was investigated in Refs. Q, Q. We will assume that expansion exists in a certain class 
of schemes. 

On the other hand, the RG- improved perturbation expansion for the observable 0(Q 2 ) is given by 



O pt (Q 2 ) = a pt (Q 2 ) + V c„_ ia "(Q 2 ) 



(4) 



n=2 



Expanding a p t(te Q 2 ) around t 



inside the integrals in Eq. must give Eq. (j^J. 



The skeleton expansion is a reorganization of the perturbation series such that each term in corresponds to 
the sum of an infinite number of Feynman diagrams. These sums, however, do not converge. Although it is possible 
to assign a value to these sums, this value is not unique, a renormalon ambiguity is present. In formulation Q the 
ambiguities arise from the (nonphysical) Landau singularities of a p t(te c Q 2 ) in the non-perturbative space-like region 
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< Q 2 < A 2 The difference between possible integration paths (prescriptions) is a measure of the size of the 
renormalon ambiguity. 

The perturbative coupling is a solution of the n-loop RG equation JSJ) . It can be found iteratively for Q 2 ^> A 2 



« P t(Q 2 ) = EE% I ^ I (5) 

i=l i=o 

where L — log(Q 2 /A 2 ) and kij are constants depending on the /3-function coefficients. At energies Q < 1 GeV the 
perturbative result 10}, (0) is not reliable. At these energies a pt starts being dominated by the Landau singularities 
at < Q 2 < A 2 . These singularities are a consequence of the perturbative RG Eq. and are located in the region 
where this equation is not valid. Furthermore, from general arguments (causality) |5| one concludes that the coupling 
parameter must be analytic in the whole Q 2 -plane excluding the time-like (Minkowskian) semiaxis. On this semiaxis, 
singularities associated with asymptotic states appear. 

With this motivation it scorns reasonable to replace the perturbative coupling by a new coupling A\{Q 2 ) differing 
from a p t(Q 2 ) significantly only in the non-perturbative region and having the required analyticity properties. This 
replacement is not unique and should be considered as a phenomcnological model. Using a dispersion relation, Shirkov 
and Solovtsov p proposed the following replacement: 

1 f°° da , _ .(MA) rr .2\ 1 I"" da 



7T J_ A 2 a + Q 2 7T J Q a + Q A 

where p\ = Im[a pt (— a — ie)], and a p t is, e.g., given by Eq. JSJl. 1 We will refer to this as the minimal analytic (MA) 
procedure. In MA the discontinuity of the analytic coupling along the Minkowskian semiaxis is by construction the 
same as the one of a p t. 2 Below we shall consider generalizations of this analytization procedure. 

Once one analytizes a p t(Q 2 ) using Eq. JJjJ or other procedure, the question arises how to treat a known truncated 
perturbation series (TPS). For perturbation series @ there is no unique way to analytize higher powers of et pt . One 
possibility is to apply the MA procedure to each power of a p t 01 : 



- 4 MA '<^» -/ ^ftW (k=l,2,...), (7) 



where pk = Im[a pt (— a — ie)} and a pt , e.g., is given by Eq. JHJ). 3 Other choices could be, e.g., a pt i-» A k 1 ,A k 1 ~ 2 A 2 , etc. 
In this paper we propose a method to analytize a TPS based on the skeleton expansion (I}, in any chosen version 
of QCD with analytic Ai(Q 2 ). In Eq. the replacement a pt i— > Ai is made, making the skeleton expansion terms 
well-defined integrals. Next, we Taylor-expand each Ai(tje c Q 2 ) there around a specific InQ 2 = \n(t r e c Q 2 ). In these 
expansions, we denote 

If we know in expansion |0J TPS with n max = 3 (i.e., c\ and C2), it is convenient to introduce the analytic couplings 
A2 and A3 according to 

A 2 m=A 2 (Ql) + ^-A 3 (Ql) , 

Po (9) 

X 3 (Q 2 ) = ^i 3 (g 2 ) . 

When replacing here Ak ^> a pt ; we obtain the corresponding truncated RG equations of perturbative QCD (pQCD). 
Thus, once one chooses a particular analytic coupling A\, the functions Ak with k > 2 are defined by Eqs. ©, or 
by a higher-n max version thereof. In MA model ©, the results 10 and © merge when n max increases (cf. also 
0). In our approach, the basic set of functions is: (Ai and its derivatives A2,As,...). The set (Ai, A2, A3, . . .) was 
introduced for the convenience of comparison with pQCD. 



1 An efficient method for evaluation of *4j MA ' was developed in Ref. 0. A different evaluation of o^t and pi was presented in Ref. 0. 

2 Other analytization procedures of a p t focus on the analyticity proper ties of the beta function [3] ■ 

3 An extension of Eq. J7J to noninteger powers was developed in Refs. Illl . 
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The afore-mentioned Taylor-expansion around InQ 2 = ln(t*e c Q 2 ) for the order n max = 3 gives 

A(fc c Q 2 ) « ^i(Q')-/3oln(tA*)X 2 (g2)+/3 2 ln 2 (t/^)X 3 (g2) (10) 
= A 1 (Ql)-poHt/t*)A 2 (Ql) + [^ln 2 (t/U)-p 1 Ht/U)]A 3 (Ql). (11) 

Keeping terms corresponding to the third-order approximation we obtain a truncated analytic version of 0{Q 2 ): 

+ [(/?oif (**) + M?{U))Az{Ql) + 2sf^ (U)A 1 (Ql)A 2 (Ql) + s° A\(Ql)\ , (12) 
where Q 2 = t*e c Q 2 and the momenta are 



f?(t*)= -F£(t)(-io g t/uy, 

Jo 1 

/£•(*•) = J ^-^-Fo(ti,t2)(-\ogt 1 /t^(-logt 2 /uy. 



(13) 



Equation l|12|) is the result of the proposed analytization procedure for the series (@J truncated at ~ a| t . In the 
perturbative region one has Ak ~ tipt and Eqs. l(T2l and merge. If _4i(<2 2 ) is well behaved at the origin then 
all Ak{Q 2 ) (k > 2) vanish at this point. This follows from the (truncated) RG-like Eqs. 10. Comparison between 
Eqs. (JT21 and 10} gives 

Cl =/3 /f(e- c ) + S f, 

We assume that the skeleton expansion coefficients s^* and characteristic functions F@(ti, . . . ,tj) are n/-independent 
when C is n/-independent. Consequently, in the class of schemes where the coefficients Cj of expansion Q are 
polynomials in u/ (<-> /3o) of order j, relations ((111) give us coefficients sf and sip and momenta /f, and /f . 
We shall consider observables for which c\ and c 2 are known. This approach can be continued to higher orders. The 
afore-mentioned class of schemes is parametrized by the RG (3j coefficients (j > 2) which are polynomials in n/ (■*-> /3o) 
of order j such that /3j = bjo + bjifto + ■ ■ ■ bjjf3y, where bjo = bjo of MS scheme and bjk (k > 1) are the free scheme 
parameters. 

However, the knowledge of the perturbation coefficients Cj by itself is not enough to obtain the higher-order co- 
efficients which are not included in Eq. 112(1 . At fourth-order, the coefficients at A 2 and ^1^3 cannot be obtained 
without certain assumptions for the characteristic function F^^i,^) (l2^ . 

The leading characteristic functio n F^ (t) is known for many observables on the basis of their all-order large-ny (<-> 
large-/?o) perturbation expansion 0, Il3|| . Therefore, we propose to keep the leading skeleton term unexpanded, but 
to expand the other terms as in Eq. 1(12(1 

{ Zl(Q 2 ) = r T f oW Mte c Q 2 ) + sfAKQl) + \s°A\{Ql) + 2 s f^ Q (t^)A 1 (Q 2 2 )A 2 (Q 2 2 )} , (15) 
Jo 1 L J 

where we used two different expansion scales for the NL and NNL skeleton terms: Q\ = t^e c Q 2 and Q 2 = t^e c Q 2 , 
respectively. Since ffo(t^) = ffo^ ) + m^ 2 -* + C and /fo^ - * 7 ) is known, it is convenient to use a scale of the 
BLM type 0, [l|: = t ( ? ] = exp(-C* - A° (e~ c )) such that f^ a (t (2) ) = 0. Consequently, the A X A 2 term in 

Eq. (|T5|l disappears. Further, the scheme dependence at this level shows up as the dependence of and on & 2 i 
and &22, respectively 01 . This allows us to fix the latter two coefficients, for each specific observable, in such a way 
that s® — and, e.g., = 1. For example, if the starting scheme is MS (621, b 22 , and C = —5/3), the new scheme 
coefficients b 2 j are 

r - 107_ 

»21 = »21 + C 2 + -TTTCll: 

19 - < 16 ) 

^22 = ^22 + c 2 i - -^cn + 2C C10, 

where Cjk are expansion coefficients of the perturbation coefficient Cj, Eq. Q, in powers of (3q, in MS scheme: 
Cj — J^oCjkPo- I n the scheme 1(16(1 . the skeleton-based expansion (11511 reduces to 

( ™\Q 2 ) = f J F$(t) A^Q 2 ) + S ?Al(e c Q 2 ) + O n , (17) 
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where O n = O4 are now formally terms of fourth order (~ Af, A 2 A2, ■ ■ ■ )• We will call formula l|17(l the first variant 
("vl") of our skeleton-motivated evaluation approach. Adopting the scheme (|16|) . or higher order generalizations of 
it, higher order contributions are absorbed in the two terms of Eq. (|17|l . This scheme- fixing method is particularly 
useful at low energies where scheme dependence is important. 

The skeleton QCD expansion, if it exists, is probably valid only in a specific (yet unknown) "skeleton" scheme 
A possible difference between the latter and the schemes used here will result in a difference in the evaluation of the 
observable 0(Q 2 ). This difference, when re-expanded in a p t, is at most ~ a^ t subleading-po (i.e., ~ f3 2 a pt ). 

The derivation up until now allows us to present yet another variant ( "v2" ) of the evaluation approach, by keeping 
the scheme (JTBJ) and simply replacing A 2 (e c Q 2 ) by Ai{e c Q' 1 ) in Eq. (|T7|) 

0^\Q 2 ) = Tt F o(t) Mte c Q 2 ) + s?A 2 (e c Q 2 ) + O n , (18) 

This formula can be obtained by repeating the previous derivation, but starting with the skeleton expansion Q 
without the analytization replacements a p t <— > Ai there. All the expansions are then obtained as previously, but with 
a™ t instead of A n , A n -iAi, etc. In this variant, the analytization is performed at the end, by replacing a p t 1— » Ai in 
the leading-skeleton integral, and replacing a 2 t i— > A2 in the term proportional to sf , leading to Eq. I|18[l . 

We wish to stress that neither variant of the evaluation approach relies on the existence of the skeleton expansion. 
Our derivation can be interpreted in the following alternative way: The formal skeleton expansion Q provides us with 
the tools to separate the perturbation series of the observable into several perturbation subseries. The first subseries 
(from the leading skeleton term) includes all the leading-/3 terms, the second subseries (from the subleading skeleton 
term) includes all the leading-/^ terms of the rest, etc. Each of these perturbation subseries is renormalization scale 
invariant. A specific renormalization scheme (/?2, p3, ■ ■ •) is then found such that all the perturbation subseries vanish, 
except the first two. In the end, the analytization of the two surviving subseries is performed. 

If the perturbation coefficient C3 is known, then the entire described procedure can be carried out to one higher 
order, i.e., the /^-coefficients b^j {j — 1, 2, 3) can be determined so that in Eqs. (fT?|l - (fT%f O n = O5 (~ Ai,A\A2, ■ • ■), 
under certain assumptions for the function f^ii,^) |T^ . For example, for the massless Adler function, the C3 
coefficient has been estimated as a polynomial in nj to a high degree of accuracy [l5| , and the scheme can be found 
such that in Eqs. (|17|) - (|18|l O n — O5. Of course, the higher order analytic couplings A n (n > 2) are defined in this 
case by the n max = 4 extension of Eqs. © 

MQl) = A 2 (Ql) + %-MQ 2 *) + ttMQH) , 
Po Po 



A 3 (Ql 

MQl) = MQ 



MQl) + l^MQl), (19) 



and expansion of Ai(te c Q 2 ) is now performed up to and including Aa(Q 2 ), in contrast to Eq. I|10JI . 

In practical evaluations, the form of the analytic coupling parameter A\{Q 2 ) has to be specified. The most straight- 
forward is the minimal analytic (MA) coupling ©. The latter model gives the value A = A( C= _ 5 / 3 ) « 0.4 GeV (in MS 
and with rif = 3) from fitting high energy QCD observables |16j . However, in order to reproduce the measured value 
of the semihadronic tau decay ratio r T , it requires introduction of heavy first generation quark masses m u « « 0.25 
GeV [13. Another possibility would be to modify the MA-coupling at low energies, e.g., in the following manner 



AM 1)fn 2,_ M 2 Q 2 M 2 Q 2 1 [°° da Pl (a)(a M 2 ) 

(Q) - Cf (Q 2 + M?) 2+k0 (Q 2 + M 2 ) + (Q 2 +M 2 )nJ a=M?t a(a + Q 2 ) ' 



In this "Ml" model, fco, c/, cq — M 2 /A 2 , and c r = M 2 /A 2 are four dimensionless and C-independent parameters 

which determine the low energy modification of the coupling (a special case, k = —I, of Ml was presented in Ref. [ljj). 

2 

In general, at high energies, this coupling differs from the MA-coupling by ~ A /Q 2 . However, requiring that the 

4 

difference be only ~ A /Q fixes parameter fco in terms of the other three. Consequently, A w 0.4 GeV from fitting to 

high energy QCD observables, as in the MA case. The remaining three parameters can be fixed by requiring that the 

experimental value of r T and some other low energy observable, e.g., Bjorken polarized sum rule dt(Q 2 ), be reproduced 

by the afore-mentioned procedure. The experimental values of these two observables are r T = 0.196 ± 0.010 |l8| and 

db(Q 2 ) — 0.16 ± 0.11 at Q 2 = 2 GeV 2 [l!j, where the normalization was chosen such that r T = a pt + 0(a 2 t ) and 

similarly for db- The quark mass effects are subtracted (not contained) here. 
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The use of the MA-coupling, in our approach (|16[) - i|17|) and with massless first three quarks and A = 0.4 GeV, gives 
db{Q 2 — 2.) « 0.13 in vl (0.14 in v2) which is acceptable, and r T — 0.140 in vl (0.139 in v2) which is not acceptable. 4 
If we require, in model Ml, the reproduction of r T « 0.196 and d b (Q 2 — 2.) » 0.13 — 0.14, in the evaluation approach 
vl, with massless quarks, and A = 0.4 GeV, we obtain for the choice Co = 2 the values of c r « 0.5 and c/ « 1.7. 
Changing Co while keeping it ~ 1 gives us by the same procedure different values of c r and Cf to reproduce the 
afore-mentioned values of r T and d b . In such cases, yet another low energy observable, the (massless) Adler function, 
remains quite stable under the variation of cq. In Ml we will take cq = 2, c r = 0.5, c/ = 1.7 (and A = 0.4 GeV). 
With these parameter values: vl approach i|17|) gives r T = 0.197 (0.202 when C„ = 4 for d{m 2 T e i9 ) in Eq. JTJl) and 
d b (Q 2 = 2.) w 0.14; v2 approach JTHJ gives r T = 0.210 (0.201 when O n = C 4 ) and <4(Q 2 = 2.) w 0.14. 

Yet another, simpler, modification of the MA-model is 

A[ M2 \Q 2 ) = 4 MA '(Q2 )+Hu ^£L- . (21) 

Q 2 + M 2 

In this "M2" model, 5 we will take the parameter values c v = 0.2, c = Mq/A 2 sa 0.56 and A = 0.4 GeV. With these 
values, vl approach (JT7J gives r T = 0.188 (0.194 when O n = £> 4 ) and c4(Q 2 = 2.) w 0.18; v2 approach {THJl gives 
r T = 0.188 (0.193 when O n = C 4 ) and d b {Q 2 = 2.) w 0.19. 

Having fixed the parameters in the afore-mentioned models Ml and M2, we present in Fig. ^ low energy results of 
the Adler function d(Q 2 ) associated with the hadronic part of the electromagnetic current, in the models MA, Ml and 
M2 (for a different evaluation of d(Q 2 ), cf. |23). The normalization was taken again such that in the massless quark 
limit for rif — 3: d(Q 2 ) = a pt + 0(a 2 t ). The lower curves in Fig. 2] represent the results of the vl-evaluation l|17|) with 
three massless quarks, in the scheme where O n = O5. The higher curves represent the full quantity, i.e., the effects of 
the massive c and b quarks are added, with the coefficients as given in Ref. |24[ (d(Q 2 ) = (1/2)Z?(Q 2 ) — 1, where D is 
defined in [24j). In the contributions of c and b, we simply replaced a p t(Q 2 ) and a 2 t (Q 2 ) by Ai(Q 2 ) and A2(Q 2 ), and 
used A = A (C = —5/3). A 2 {Q 2 ) was constructed by the fourth-order relations JT5J|. The indicated ± uncertainties 
in these curves are those charm contributions which are oc A 2 (Q 2 )- In contrast to the massless contributions, we do 
not have yet a systematic way to analytize the massive quark contributions. The experimental results [24[ and the 
truncated pQCD series were included for comparison. Figure 1 shows that analytic versions of QCD (MA, Ml, M2) 
in conjunction with the skeleton-motivated approach l|17|l give results that at low energies Q ~ 1 GeV behave much 
better than pQCD. 

Application of variant 2 of our approach, Eq. (|18f) . gives for the Adler function results which are very close to those 
of variant 1, Eq. (|17fl . because the coefficient sf is small: sf = 1/12. For the Bjorken polarized sum rule db(Q 2 ), 
this is not so, because sf = —11/12 is appreciable. Therefore, we will take the Bjorken sum rule as a case to look at 
numerical differences between evaluations in our approaches vl and v2, Eqs. H17fl - I|18|) . I n addition, we will compare 
with the approach of Milton et al. and Shirkov [TflEl (MSSSh) 

OmsssUQ 2 ) = A^Q 2 ) + Cl A 2 (Q 2 ) + c 2 A 3 (Q 2 ) . (22) 

In principle, the comparison between the approaches can be carried out in any anaylitic QCD model for the coupling 
parameter. However, the MSSSh approach has been applied in the literature in the MA model. Therefore, we carry 
out the comparison in this model. In Fig.|21we present the MA-model predictions for db(Q 2 ), with various evaluations, 
at third order [O n = O4 in Eqs. (|T7Jl - ljT%|l ]. The results of the MSSSh approach are presented in two schemes: MS 
("bMS"); and in the "b2Sk" scheme where 02 is determined by Eqs. ltTH|l for d b . The MSSSh approach uses the 
three-loop MA-expressions ^(Q 2 ) and As(Q 2 ) of Eq. Q (Refs. |10L Il6| ). and variant 2 of our approach, Eq. i|18fl . 
uses A 2 {Q 2 ) from Eqs. ©. Figure [3 shows that the evaluation of d b {Q 2 ) with variant 1 of our approach Eq. I|17|) 
("Sk vl b2Sk") gives, at low energies, results which differ significantly from the MSSSh approach. On the other hand, 
variant 2 of our approach ("Sk v2 b2Sk"), i.e., Eq. (|18() . gives results which are, apparently accidentally, very close to 
those of the MSSSh approach in the MS scheme. 

In summary, we presented two variants of skeleton-expansion- motivated evaluation of observables in analytic versions 
of QCD, Eqs. (|17fl - l|18fl . The first variant follows more closely the skeleton expansion, in the sense that the analytization 



4 For the corresponding massless Adler function d(Q 2 ) we use the scheme where in Eqs. {T^-Eg) O n = O5 E3- 

If taking a scheme with 

by one order lower precision (O n = O4, /3s 1— > 0), the value changes to r T = 0.146 in vl (0.145 in v2). Observable r T is evaluated by 
first evaluating the Adler function d(Q 2 ) for complex values Q 2 = m 2 e lS and then applying the standard contour integration in the 
Q 2 -plane. The function Fq (t) for d(Q' 2 ) was obtained in Ref. Q. F$(t) for d b (Q 2 ) can be obtained from the known large-ra^ expansion 
of d b (Q 2 ) |2(1 , using the technique of Ref. , and the full perturbation coefficients c\ and C2 from |21| . 

5 In Ref. |22| . power correction terms 1/(Q 2 )" were added to ,4^ MA '(Q 2 ), but with a somewhat different motivation. 
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is performed at the beginning, in the skeleton expansion (a pt i— > Ax). In the second variant, the analytization is 
performed at the end, in the form a^ t i— > .At. The second variant can be regarded as a generalization of the evaluation 
approach of Milton et al. and Shirkov llfil] (MSSSh), now including the leading- (3$ terms to all orders in the 
coupling parameter. Both variants use the formal structure of the skeleton expansion in order to divide the original 
perturbation expansion into a sum of subseries (skeleton terms), each of them renormalization scale invariant, and 
then using a scheme where only the first two subseries survive. Further, we introduced two alternative models (Ml, 
M2) of analytic QCD for the coupling parameter which, for certain values of the model parameters, reproduce the 
measured values of the semihadronic r decay ratio. 
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FIG. 1: Adler function as predicted by pQCD, and by our approach in several analytic QCD models (see the text). 
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FIG. 2: Bjorken polarized sum rule in MA as predicted by two variants of the approach lllllKl (MSSSh), and by our approach. 



